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ABSTRACT 


The  Naval  Postgraduate  School  simulation  model,  QDYN85, 
was  used  to  investigate  the  sputtering  of  nitrogen  reacted 
tungsten  and  molybdenum  targets,  where  the  nitrogen  was 
placed  in  equilibrium  positions  both  above  and  below  the 
target  surface  plane.  The  cascades  were  initiated  by 
bombarding  the  targets  with  argon  ions  with  energies  varying 
from  1  kev  to  3  kev.  The  simulation  results,  when  compared 
to  experimental  results  obtained  by  Winters,  indicate  that 
the  nitrogen  position  is  not  below  the  surface  plane  for 
either  target  but  is  located  near  the  point:  where  the  adatom 
would  be  equidistant  from  its  five  nearest  neighbors. 
Furthermore,  the  results  show  no  multimer  formation  as  a 
result  of  lattice  fragmentation  but  do  show  dimer  formation 
as  a  result  of  recombination  above  the  target  of 
individually  sputtered  atoms. 
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I.  OVERVIEW 


A.  DEFINITION  AND  EARLY  EXPERIMENTAL  WORK 

Sputtering  is  the  ejection  of  surface  atoms  from  an 
object  due  to  energy  transfer  from  an  energetic  projectile 
ion  to  the  atoms  of  the  object.  During  each  collision 
between  the  energetic  ion  and  a  target  atom,  some  energy  and 
momentum  is  given  up  by  the  ion.  If  this  transferred  energy 
exceeds  the  atoms’  lattice  binding  energy,  the  struck  atom 
(recoil  atom)  will  be  displaced  out  of  its  lattice  site  and 
will  collide  with  other  target  atoms  initiating  a  cascade. 
If  a  recodl  atom  is  located  at  the  surface  of  the  target, 
and  its  acquired  energy  exceeds  the  surface  binding  energy, 
and  the  direction  of  motion  is  away  from  the  target,  that 
atom  will  be  ejected.  This  cascade  sequence  continues  until 
the  initial  energy  has  been  dissipated  so  far  that  further 
collisions  could  not  possibly  result  in  target  atom 
displacements . 

The  sputtering  process  was  first  reported  in  1852  by 
GROVE  (Ref  1)  when  he  noticed  a  deterioration  of  cathodes 
and  a  blackening  of  the  glass  of  glow  discharge  tubes.  Over 
the  next  several  decades,  a  considerable  amount  of  work  was 
done  in  the  field  of  sputtering.  However,  a  lack  of 
understanding  of  how  various  parameters  such  as  ambient 
pressure,  target  purity  and  angle  of  incidence  of  the  ion 
beam  would  affect  the  sputtering  yield  ,Y,  (  defined  as  the 


number  of  ejected  particles  per  incident  ion  )  resulted  in 
non-reproducible  experiments.  This  effectively  negated  these 
early  efforts. 

By  1902,  GOLDSTEIN  (Ref  2)  presented  evidence  that 
sputtering  was  caused  by  positive  ions  impacting  on  the 
metal  cathode.  Since  this  was  considered  a  deterioration  and 
thus  an  undesireable  effect  at  the  time,  work  in  the  area 
was  centered  on  controlling  the  sputtering  effect.  However, 
some  early  researchers  such  as  PLUCKER  (Ref  3),  recognized 
the  possibilities  for  other  experimental  work  of  the  effect 
in  thin  film  manufacturing  and  removal.  PLUCKER  and  other 
pioneers  in  sputtering  laid  the  foundation  which  established 
sputtering  as  a  field  of  quantitative  experimental  research. 

B.  THEORETICAL  AND  EXPERIMENTAL  DEVELOPMENTS 

The  dependence  of  the  sputtering  yield  on  the  ambient 
pressure  was  first  demonstrated  by  PENNING  and  MOUBIS  (Ref 
4)  in  1940.  They  realized  that  the  number  of  collisions 
between  atoms  ejected  from  the  target  and  atoms  of  the 
ambient  atmosphere  of  the  vacuum  system  increased  when  the 
pressure  was  increased.  They  showed  that  this  was  causing 
backscattering  of  sputtered  atoms  into  the  target .  To 
overcome  this  effect,  they  kept  the  ambient  pressure  no 
higher  than  10-5  Torr.  At  this  level,  reproducible  results 
were  obtained  for  ion  energies  greater  than  500  ev.  For  ion 
energies  less  than  500  ev,  the  sputtering  action  was 
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insufficient  to  maintain  a  'clean"  surface. 

STARK  (Ref  5  ),  pursuing  sputtering  as  an  individual, 
atomic  scale  event,  developed  two  separate  models  of  the 
process.  His  hot  spot  model  was  based  on  the  assumption  that 
the  evaporation  of  target  material  from  a  very  small  area 
was  caused  by  the  high  temperature  created  by  the  ion’s 
impact.  VON  HIPPLE  (Ref  6)  further  developed  this  approach 
and  expended  a  considerable  amount  of  effort  attempting  to 
formulate  a  sputtering  theory  based  on  this  premise,  stating 
that  this  approach  was  the  only  feasible  avenue  to  address 
the  complex  statistics  of  the  collision  processes  occuring 
in  a  sputtering  event.  It  should  be  noted  that  although  this 
approach  has  recently  been  reexamined  by  KELLY  and  coworkers 
(Ref  7),  it  has  not  enjoyed  the  attention  and  development 
accorded  to  Stark's  other  model,  the  "collision"  model. 
In  1921,  KINGDON  and  LANGMUIR  (Ref  8)  produced  a  successful 
application  of  Stark’s  collision  theory  to  the  analysis  of 
ion  induced  desorption  of  monolayers. 

In  the  1950’ s,  KEYWELL  (Ref  9,10)  attempted  to  formulate 
Stark's  collision  model  in  terms  of  neutron  transport 
theory.  This  and  subsequent  calculations  by  HARRISON  (Ref 
11)  were  significant  because  probability  concepts,  in  the 
form  of  collision  cross  sections,  were  introduced  into  the 
theory.  Although  important  steps,  these  theories  were 
not  useful  because  they  depended  on  a  large  number  of 
unknown  parameters;  so  quantitative  yield  values  could  not 


be  calculated.  In  1956,  WEHNER  (Ref  12)  demonstrated  the 
effect  that  the  target1 s  crystalline  structure  could  have  on 
the  sputtering  process  when  he  observed  spot  patterns  while 
bombarding  various  monocrystals .  These  spots  were  atom 
ejection  patterns  which  corresponded  to  the  close-packed 
directions  of  the  monocrystal  target. 

Futher  theoretical  work  by  SILSBEE  (Ref  13)  led  to  the 
suggestion  that  energy  was  transported  through  a  lattice 
structure  by  momentum  transfer  along  close  packed 
directions.  This  phenomenon  is  called  focusing  and  the 
transport  mechanism  is  called  a  f ocuson.  Extensions  of 
Silsbee’s  focusing  concept  were  developed  by  THOMPSON  (Ref 

14)  and  by  experimental  work  by  NELSON  and  THOMPSON  (Ref 

15) .  While  the  focuson  mechanism  exists,  subsequent 
investigations  have  shown  that,  its  contribution  to 
sputtering  is  small  and  that,  by  itself,  it  could  not 
explain  the  spot  patterns  observed  by  WEHNER.  By  analogy, 
ROBINSON  and  OEN  (Ref  16)  expected  that  a  channelling  effect 
on  an  ion  through  the  lattice  would  be  found.  This 
channelling  effect  would,  of  course,  be  dependent  on  the 
angle  of  incidence  and  the  crystalline  structure.  They 
effectivly  discovered  this  while  studying  ion  penetration  of 
a  crystalline  structure  using  a  simulation  model.  This 
phenomenon  was  confirmed  later  experimentally  by  DAVIES  et 
al.  (Ref  17) 

Alternative  theoretical  approaches  to  the  sputtering 


problem  culminated  in  the  theories  of  SIGMUND  (Ref  IS)  and 
THOMPSON  (Ref  14)  which  consider  ensemble  averages  over 
distribution  functions.  This  approach  suppresses  any  use  of 
a  target  lattice  structure  in  favor  of  an  amorphous  target. 
At  higher  energies,  reasonably  correct  values  of  the 
sputtering  yield,  the  correct  form  for  the  ejected  atom 
energy  distribution  function  and  the  ejected  atom  angular 
distribution  function  from  polycrystalline  targets  are 
obtained.  This  theory  is  relativly  inaccurate  at  low  ion 
energies  since  it  is  in  this  region  that  lattice  effects  are 
most  pronounced  and  are  of  significant  consequence.  While 
refinements  and  improvements  followed,  attempts  to  reduce  a 
many  -  bodied,  multiple  interaction  event  into  an 
analytically  feasible  model  have  still  not  been  able  to 
fully  explain  observed  behavior. 

C.  COMPUTER  SIMULATION 

With  the  advent  of  high  speed,  sophisticated  computers, 
it  became  possible  to  simulate  a  sequence  of  atomic 
collision  events.  The  development  of  the  simulation  tool 
proceeded  along  two  paths.  Since  the  development  of 
sputtering  as  a  series  of  independent  binary  collisions  was 
a  natural  outgrowth  of  STARK’S  initial  theories,  this  type 
of  simulation  was  developed  first.  One  of  the  principle 
philosophical  drawbacks  to  this  approach  is  that  a  complete 
understanding  of  the  physical  model  is  assumed  and  that  the 


assumption  that  each  event  is  separate  from  another  can  be 
justified.  While  this  concept  may  be  practical  for  higher 
initial  energies,  it  is  suspect  for  lower  energies  where 
multiple  collisions  are  important. 

Starting  in  1960,  a  different  approach  to  the  sputtering 
simulation  developed.  It  was  recognized  that,  particularly 
at  low  energies,  a  cascade  could  not  be  treated  as  a  series 
of  independent  events.  Instead,  a  multiple  interaction  (MI) 
procedure  was  required.  In  this  type  of  procedure,  during 
a  specific  short  period  of  time,  many  atoms  of  the  target 
will  simultaneously  affect  the  cascade's  development.  This 
type  of  simulation  was  introduced  by  GIBSON,  GOLANB,  MILGRAM 
and  VINEYARD  (Ref  19)  when  they  developed  a  computer 
simulation  model  of  the  motions  of'  copper  atoms  in -a 
crystalline  target  after  one  of  the  atoms  had  been  struck  by 
an  incident  ion. 

This  time-step  approach  to  sputtering  simulation  has, 
over  the  years,  illuminated  important  aspects  of  the  process 
which  were  masked  by  the  ensemble  averaging  theories.  The 
realizations  that  sputtering  occurs  primarily  in  the  first 
three  layers  of  the  target  and  that  focusons  do  not,  to  a 
great  extent,  affect  the  sputtering  yield  were  made  in  this 
manner.  This  thesis  will  use  the  MI  program  QDYN86  to 
simulate  the  problem  outlined  below. 

While  a  computer  simulation  can  provide  insight  into  the 
sputtering  process,  it  is  just  that;  a  simulation.  An 
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understanding  of  the  physical  situation  must  be  attained  and 
balanced  by  physics  and  mathematics  before  a  simulation  can 
be  successfully  achieved.  Some  fundamentals,  such  as  the 
interatomic  potential  functions,  are  still  not  fully 
understood.  In  this  respect,  computer  simulations  cannot  and 
will  not  supercede  experimentation.  They  can  however, 
provide  support  and  explanation  of  observed  phenomenan  and, 
more  importantly,  be  used  as  a  cutting  edge  of  science  in 
the  probing  of  new  ideas  and  concepts. 

D.  APPLICATIONS 

Sputtering  simulations  are  not  limited  to  monocrystalline 
targets.  Indeed,  much  of  the  effort  currently  being  expended 
lies  in  the  arena  of  alloys  and  chemically  reacted  surface 
bombardment.  Chemical  reaction  of  the  surfaces  involves  the 
introduction  of  an  additional,  different  atom  onto  the 
surface  of  a  crystalline  target.  This  additional  atom  is 
called  an  adatom .  The  locations  of  these  adatoms,  prior  to 
bombardment,  is  usually  unknown,  however,  for  various 
crystalline  structures,  anticipated  locations  may  be 
described.  Both  tungsten  and  molybdenum  are  BCC  lattice 
structures  and  the  four-fold  position  (Fig  1 )  is  a  possible 
candidate,  depending  on  the  amount  of  coverage  desired.  In 
that  location,  the  adatom  would  sit  in  a  position  between 
the  four  face  atoms  and  a  center  atom  such  that  either  the 
forces  acting  on  it,  such  as  lattice  and  surface  binding 
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energies,  equalize,  or  it  could  be  in  a  position  such  that 
the  distance  to  it’s  nearest  neighbors  would  be  equal  (Fig 
1).  Whether  the  adatom  is  more  likly  to  be  located  above  or 
below  the  surface  plane  and  which  of  the  above  situations  is 
applicable  is  an  unanswered  question  for  most  systems. 

The  realization  that  important  real  world  applications 
can  be  addressed  in  simulations  has  spurred  their 
application  to  surface  science  investigations.  Corrosion, 
the  formation  of  oxides  on  the  surface  of  a  substance,  can 
be  studied.  Possible  solutions  to  its  destructive  effects 
such  as  adatom  placement  to  either  prohibit  or  slow  down  the 
introduction  of  the  oxygen  atom  into  the  lattice  can  be 
simulated.  Conversly,  positioning  may  allow  the  adatom  to 
act  as  a  catalyst,  enhancing  a  desired  reaction. 


* 
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II .OBJECTIVES 


Winters  (Ref  20)  has  conducted  experiments  investigating 
sputtering  cross  sections  of  both  tungsten  and  molybdenum 
targets  which  have  been  reacted  with  nitrogen  adatoms,  when 
bombarded  by  various  ions.  These  experiments  showed  that  the 
nitrogen  sputtering  yield  tended  to  increase  as  the  atomic 
weight  of  the  target  atoms  increased.  Thus,  assuming  a  (001) 
face  on  the  target,  the  sputtering  yield  of  the  nitrogen 
adatoms  was  reasoned  to  be  dominated  by  the  mass  ratio  of 
the  substrate  to  the  adatoms.  Meyerhof f  (Ref  21),  using  a 
simulation  model,  determined  that  the  distance  between  the 
absorbate  and  the  substrate  atoms  was  more  significant  than 
the  mass  ratio. 

Both  Winters  and  Meyerhoff  assumed  that  the  location  of 
the  nitrogen  adatoms  was  in  the  four-fold  position  (Fig  1). 
Meyerhoff  conducted  his  simulation  assuming  a  nitrogen 
location  slightly  above  the  suface  plane,  resulting  in  a 
position  of  equal  distance  between  the  adatom  and  all  five 
nearest-neighbors .  Earlier  work  by  Calvenna  and  Schmidt  (Ref 
22)  and  Adams  and  Germer  (Ref  23)  also  assumed  the  four-fold 
position  however,  there  was  little  regard  •  given  to  target 
atom  displacement  due  to  the  introduction  of  the  nitrogen 
adatoms.  Griffiths  et .  al.  (Ref  24)  incorporated  this  idea 
in  a  series  of  experimants  using  low  energy  electron 
diffraction  (LEED).  Using  the  observed  patterns  for  various 


coverages  and  primary  beam  energies,  they  were  able  to  show 
that,  at  a  fractional  coverage  (  defined  as  the  ratio,  for  a 
specific  face  of  the  target,  of  the  number  of  adatoms  to  the 
number  of  substrate  atoms.)  of  0.5  and  a  surface  temperature 
of  1000° K,  the  nitrogen  adatoms  were  unstable  and  were 
rapidly  absorbed  into  the  bulk.  If  the  coverage  was  reduced 
to  approximatly  0.4,  they  found  that  temperatures  up  to 
HOOoK  did  not  result  in  any  further  absorption.  A  coverage 
of  0.4  however,  required  the  formation  of  "islands"  of 
substrate  with  nitrogen  adatoms  on  the  target  (Fig  2).  The 
target  used  in  this  thesis  can  be  interpreted  as  one  of 
these  "islands".  To  investigate  this  type  of  structure, 
various  laser  diffraction  gratings  were  constructed  and 
tested  by  Griffiths  et  al .  A  grating  with  a  random  number  of 
nitrogen  adatoms,  but  with  an  average  island  size  of  4X4, 
produced  a  pattern  similar  to  that  observed  when  LEED  beams 
were  used  on  an  actual  target  although  the  intensities  of 
various  parts  of  the  pattern  were  not  similar.  Harrison  (Ref 
25)  suggested  that  the  nitrogen  adatoms  could  be  located 
below  the  surface  close  to  an  interstitial  site.  As  stated 
earlier,  Meyerhof f  conducted  his  simulation  assuming  the 
equilibrium  position  of  the  nitrogen  adatom  above  the 
surface  plane  of  the  substrate.  This  thesis  will  study 
equilibrium  points  both  above  and  below  the  substrate 
surface  plane. 

The  primary  objective  of  this  thesis  is  to  study 


sputtering  from  BCC  crystalline  targets  of  molybdenum  and 
tungsten  which  have  had  nitrogen  adatoms  placed  in  different 
equilibrium  positions,  above  and  below  the  substrate  surface 
plane  for  both  molybdenum  and  tungsten  targets.  The 
incident  ions  will  be  argon  of  various  energies. 


1 
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III.  Simulation  Model  Development 

A.  QDYN85 

The  simulation  program  used  for  this  thesis  research 
was  QDYN85.  This  program  is  an  MI  simulation  using  a  time 
step  approach.  Initial  inputs  of  crystalline  structure, 
atomic  masses,  potential  function  parameters,  adatom 
locations,  ion  energies  and  impact  point  are  required.  Using 
Hamilton’s  equations  of  motion,  the  program  initiates  and 
develops  the  subsequent  cascade  through  the  crystal  until 
the  energy  of  the  atoms  is  lowered  to  the  point  where  no 
futher  ejections  can  occur. 

After  the  initial  ion  impact,  position,  velocity 
directions  and  energies  of  affected  target  atoms  and  adatoms 
are  computed  every  timestep.  To  compute  these  values  for 
every  atom,  every  timestep,  the  actual  program  running  time 
(cost)  would  be  excessive.  Therefore,  forces  are  not 
computed  on  a  particular  target  atom  until  it  is  actually 
hit  by  a  moving  atom.  Additionally,  the  program  maintains  a 
listing  which  is  updated  periodically,  of  each  atom’s 
nearest  neighbors . 

For  this  particular  simulation,  the  timestep  is  variable 

and  is  determined  by  a  specified  distance  divided  by  the 

highest  atomic  velocity  in  the  target.  The  distance  choosen 

for  this  calculation  is  0.1  lattice  units  (LU).  A  lattice 
* 

unit  is  defined  as  one  half  the  lattice  parameter  of  any 


cubic  crystalline  structure.  The  lattice  unit  was  created  to 
facilitate  calculations  concerning  crystal  structure  by 
eliminating  the  need  to  input  a  large  number  of  new 
variables  when  there  is  a  change  in  the  elements  being 
studied.  To  determine  the  appropriate  timestep  more 
accurately,  the  rough  "distance  over  velocity"  calculation 
is  modified  by  limiting  factors  which  take  into  account  the 
previous  velocity  thus  ensuring  a  smoother  transition  in 
time . 

Harrison  and  Jakas  (Ref  26)  present  a  more  detailed 
description  and  analysis  of  the  MI  code  of  QDYN85,  as  well 
as  a  comparison  of  its  results  with  those  obtained  by 
TRIMSP,  a  binary  collision  simulation. 

B.  POTENTIAL  FUNCTIONS 

Of  critical  importance  to  the  sputtering  process  is  the 
modeling  of  how  the  atoms  of  the  crystal  interact  with  each 
other.  To  date,  no  single  potential  function  has  been 
developed  that  fully  and  accuratly  reflects  every  sputtering 
experiment.  While  newer  potential  functions  may  be  developed 
or  the  parameters  changed  to  fit  the  observed  results,  the 
basic  mechanisms  remain  stable.  Four  basic  potential 
functions  are  used  here. 

1 . Born-Mayer 

This  is  a  strictly  repulsive  function  which  is 
considered  adequate  for  intermediate  separations.  It  has  the 
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form: 


V( r ) =  a  exp  ( -br ) . 

2 .  Moliere 

This  is  also  a  repulsive  function  but  of  the  form: 

V( r )  =  ( Zi  Zz e2 /r ) [0 . 35exp( -0 . 3r/a ) +0 . 55exp( -1 . 2r/a ) 

+0 . lexp( -6 . Or /a) ] , 

where  "a"  is  the  Firsov  length  and  is  calculated  as  follows 

a=0.8853ab/(Zii/2+Z2i/2>  )2/3, 

where  ab  is  the  Bohr  radius.  Because  the  Moliere  function 
was  originally  intended  as  an  approximation  to  the 
Thomas-Fermi  function.it  is  accepted  that  the  value  obtained 
for  "a"  may  be  modified.  This  initiated  a  new  terminology, 
the  "modified  Moliere"  potential  function. 

3 .  Morse 

The  Morse  is  a  potential  function  which  is  both 
attractive  and  repulsive  depending  on  the  separation 
distance  "r".  It  has  the  following  form 

V(r)=Deexp[-26c(r-re  )]  -  2Deexp[-ft(r-re  )  ]  , 

where  De  is  the  well  depth,  re  is  the  equilibrium  separation 
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of  an  atomic  pair  and  alpha  is  a  parameter  controlling  the 
shape  of  the  well 
4 • Morse-Moliere 

This  function  is  a  composite  potential  function 
consisting  of  sections  of  both  the  Morse  and  the  Moliere 
functions.  The  two  are  joined  together  at  a  specific  section 
by  a  cubic  spline.  The  intersection  is  determined  by  varying 
the  alpha  and  "a"  parameters  of  the  Morse  and  Moliere 
functions  respectivly  to  obtain  an  intersection  where  the 
slopes  are  as  close  together  as  possible.  The  result  is  a 
smooth  joining  with  the  repulsive  wall  governing  the 
collision  dynamics  and  an  attractive  well  which  governs  the 
sputtering. 

Obviously,  values  obtained  by  these  functions  are 
dependent  on  the  elements  used  for  the  target,  adatom  and 
the  incident  ion.  Table  1  gives  the  pertinent  data  for  the 
elements  used  in  this  simulation.  The  value  ao  is  the 
lattice  parameter,  B  is  the  binding  energy  of  nitrogen  to 
the  substrate  and  C  is  the  cohesive  energy. 

Table  Is  Natural  Data  for  Elements  Used. 


At .  Mass 

Structure 

ao  (A) 

B(eV) 

C(eV) 

Mo 

95.94 

BCC 

3.147 

N/A 

6.820 

W 

183.85 

BCC 

3.165 

N/A 

8.900 

Ar 

18.0 

N/A 

N/A 

N/A 

N/A 

N 

14.0 

N/A 

N/A 

6.5 

N/A 

5 .  Substrate-Substrate  Function 

The  potential  function  used  for  the  crystalline 
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structure  is  a  composite  Morse-Moliere  tunction.  Th^ 
development  of  the  potential  parameters  for  this  interaction 
will  be  essentially  the  same  for  the  substrate-adatom, 
substrate-ion  and  adatom-ion  interactions  described  below. 

As  stated  earlier,  the  alpha  and  the  "a"  values  were 
choosen  to  match  the  respective  potential  graph  slopes.  The 
value  of  De  was  choosen  so  that  the  energy  from  the  atoms  of 
the  lattice,  measured  at  the  center  of  the  lattice,  equaled 
the  cohesive  energy  listed  in  table  1.  The  value  of  Re  is 
the  nearest  neighbor  separation  distance.  The  values  of  Ra 
and  Rb , the  spline  boundaries,  were  determined  by  the  final 
graphs  of  the  potentials  which,  in  turn,  were  dictated  by 
the  alpha  and  "a"  values.  There  is  also  a  parameter  Rc  which 
is  used  in  the  program.  It’s  function  is  to  limit  the 
distance,  in  lattice  units,  over  which,  the  force  and 
potential  calculations  are  made.  This  value  is  1.71  LU 
except  for  the  substrate-substrate  interaction  in  which  case 
Rc  is  2.40  LU. 

6.  Adatom-Substrate  Function 


The  parameters  for  this  function  are  determined  in 
the  same  manner  as  the  substrate-substrate  function  with  the 
exception  of  the  Ra  and  De  values.  In  this  instance,  the 
value  Re  is  the  measured  separation  distance  vertically 
between  the  cubic  center  and  the  adatom  and  De  is  the 
binding  energy  of  the  nitrogen  atoms  to  the  surface, 
6 . 5eV/atom. 


7 .  Ion-Substrate/Adatom  Functions 

In  addition  to  these  adatom-substrate  potential 
functions,  others  must  be  considered.  Since  nitrogen  has 
been  introduced  into  the  target,  the  possibility  exists 
that,  at  some  time,  two  nitrogen  atoms  may  become  neighbors 
or  even  collide.  The  potential  function  used  to  describe 
this  evolution  is  a  pure  Morse  function  with  the  following 
experimentally  determined  values: 

6(.  =  2.700  i  , 

Re=  1.098  j|, 

De=  7.373  eV, 

During  the  simulation  of  a  cascade,  a  single  incident  ion 
is  injected  onto  the  target.  Since  it  is  unlikly  that  the 
argon  ion  will  react  with  the  nitrogen,  molybdenum  or  the 
tungsten,  these  potential  functions  are  assigned  as 
"modified"  Moliere’s. 

The  overall  initial  governing  factors  in  the 
determination  of  the  above  values  in  this  and  the  two 
previous  sections  were  the  desire  to  obtain  a  smooth  force 
curve  and  for  the  minimum  point  of  the  well  depth  to  occur 
when  the  adatom  was  below  the  surface  plane.  Tables  2  and  3 
show  the  comparison  between  the  values  assuming  an  adatom 

f 

{ 

location  of  0.05  LU  below  the  crystalline  plane,  and 
Meyerhoff’s.  Table  4  shows  the  values  assuming  various  above 
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plane  positions  but  incorperating  the  interatomic  forces. 


Table  2:  POTENTIAL  FUNCTION  PARAMETERS. 


(PREVIOUS) 

De(eV)  Re  ( A ) 

aiA-1 ) 

MOD  a 

Ra  (LU) 

Rb ( LU ) 

N-Mo 

1.612  1.950 

2.200 

0.0 

.910 

.  970 

N-W 

1.300  2.380 

1.650 

0.0 

1 .140 

1  .  140 

Mo-Mo 

.997  2.800 

1 . 500 

0.0 

.790 

.  830 

W-W 

1.335  2.049 

1.200 

0.0 

1 . 830 

1.900 

Rb  ( LU ) 

.  560 
.  670 
.  830 
1.130 


Rb 

(LU) 

.  750 
.  520 
.  850 

graphs 

for  the  new  values  of  the  crystal  potential  function  as 
compared  to  those  used  previously  by  Meyerhof f.  A  review  of 
the  force  graphs  using  Meyerhoff's  data,  shows  a  possible 
discontinuity  point  in  the  area  where  the  splining  between 
the  Morse  and  Moliere  potentials  occurs.  It  should  be  noted 
that,  at  the  time  of  the  research,  the  effects  the 
interatomic  forces  generated  would  have  on  the  final  results 
was  not  recognized.  The  differences  are  not  significant. 
Figures  11  through  14  show  the  plots  for  the  old  and  new 


Table  3:  POTENTIAL  FUNCTION  PARAMETERS  (NEW) 


(BELOW  PLANE) 
De(eV)  Re  (A) 

dWA’1  ) 

MOD  a 

Ra  (LU) 

N-Mo 

2.970 

1 . 524 

2 . 600 

-.800 

.  540 

N-W 

3.292 

1.524 

2.850 

-.800 

.650 

Mo -Mo 

.  997 

2.800 

1 . 519 

0.0 

.  790 

W-W 

1.335 

2.894 

1.200 

0.0 

1.120 

TABLE  4:  CRYSTAL  POTENTIAL  FUNCTION  PARAMETERS 
ABOVE  PLANE  LOCATION 


LOCATION 

De 

Re 

o< 

MOD  a 

Ra 

(eV) 

(A) 

(A*1  ) 

(LU) 

Mo-N 

.380 

1 . 600 

i  .950 

2.150 

0.0 

.  720 

Mo-N 

.  245 

1 . 997 

1.738 

2.143 

-.8 

.  500 

W-N 

.  487 

1 .712 

1.960 

2.320 

.9 

.  830 

Figures 

3-10 

show 

the  potential 

and 

force 
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Mo -Mo  and  W-W  forces.  Figure  15  shows  the  graph  of  the 
nitrogen-nitrogen  potential  function  and  figure  16  shows  the 
graphs  of  the  ion-adatom  and  the  ion-substrate  potential 
functions . 

8 . Vacuum  Phase  Potential  Parameters 

The  probability  exists  that  interactions  above  the 
surface  of  the  target  will  occur.  The  potentials  used  for 
these  interactions  are  all  pure  Morse.  There  is  a  paucity  of 
data  for  the  Mo-Mo,  W-W,  Mo-N  and  W-N  systems  and,  as  such, 
approximations  had  to  be  made.  The  values  required  are  De  , 
Re  ,  "omega1'  and  "omega-x"  where  the  later  two  are  functions 
of  the  vibrational  and  rotational  energy  of  the  system.  The 
value  <54  is  also  required  and  is  calculated,  after 
determining  De ,  as  follows: 

<54  =  ( 1 . 361 16x10“  3  )(  omega  )( m/De  ) 1 /2 

m  =  reduced  mass  of  the  system 

No  data  could  be  found  for  the  Mo-Mo,  W-W,  Mo-N  or  the  W-N 
systems.  Therefore,  as  a  compromise,  using  reference  27, 
systems  were  found  that  were  as  close  to  the  investigated 
system  as  possible  in  reduced  mass  and  in  location  on  the 
periodic  table.  For  Mo-Mo,  Ag-Br  was  used  and  VO  was 
substituted  for  Mo-N.  Similar  substitutions  were  used  for 
the  tungsten  system.  The  N-N  system  was  available.  Table  5 


summarizes  the  values  used. 


TABLE  5:  VACUUM  PHASE  POTENTIAL  PARAMETERS. 


De 

Re 

omega 

omega-x 

Mo-Mo 

1.700 

1.790 

2.390 

247.700 

.  679 

N-N 

5.760 

3.359 

1.098 

2358.570 

14.320 

Mo-N 

6.100 

1 . 965 

1 . 524 

1020.000 

4.700 

W-W 

3.100 

1.450 

2.470 

308.399 

.  960 

W-N 

4.900 

2.159 

1.738 

967.000 

4.850 

C.  TARGET  and  IMPACT  AREAS 

The  size  of  the  target  is  determined  by  the  energy  of  the 
incident  ion.  If  the  selected  target  size  is  insufficient, 
it  is  possible  that,  during  the  course  of  a  cascade,  atoms 
which  have  escaped  from  the  sides  or  the  bottom  of  the 
target  would  have  caused  additional  sputtering  to  occur  if 
the  target  had  been  larger.  Such  an  event  is  called  a 
"failure  of  containment" .  Although  a  large  target  is 
desirable  to  maximize  the  sputtering  yield,  there  exists  a 
practcal  upper  bound.  For  a  given  incident  energy, 
increasing  the  size  of  the  target  beyond  some  point  would 
waste  computer  resources.  Therefore,  a  balance  must  be 
struck  between  the  two  requirements.  Using  a  copper  target 
bombarded  with  argon,  it  was  found  that  a  target  size  of 
23X8X23  for  an  incident  energy  of  2  keV  would  provide  almost 
complete  containment. 

For  this  size  target,  the  introduction  of  60  nitrogen 
adatoms,  each  in  the  four-fold  position,  would  result  in  a 
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fractional  coverage  of  .417.  Figure  17  shows  the  (001)  face 
for  such  a  target.  Although  nitrogen  is  the  adatom  used  in 
this  thesis,  other  elements  may  be  used.  This  could  result 
in  the  adatoms  occupying  a  different  location  on  the  target. 
These  locations  may  be  the  two-fold  (equidistant  between  two 
nearest-neighbor  substrate  atoms)  or  the  A-top  (directly 
above  a  substrate  atom)  positions  (figure  18). 

A  total  of  600  impact  points  on  the  surface  were  examined 
for  each  incident  ion  energy  value.  These  impact  points,  as 
measured  on  the  x  and  z  axis,  were  varied  about  a  central, 
pre-selected  point  on  the  target.  There  were  two 
pre-selected  points,  (11.000,  11.000)  and  (13.000,  15.000) 
(Fig  17).  The  first  location  does  not  have  an  occupying 
adatom  while  the  second  does.  This  method  results  in  a  large 
number  of  independent  runs  for  a  specific  ion  energy  and 
therefore,  more  accuracy  in  the  results. 


IV.  RESULTS  AND  DISCUSSIONS 


A.  ANALYZING  AND  PLOTTING  PROGRAMS 

ANMOL  is  the  analysing  program  currently  being  used  for 
sputtering  research  at  the  Naval  Postgraduate  School.  The 
results  obtained  from  QDYN85  are  used  as  input  with  the 
ANMOL  program  formatting  and  tabulating  the  results.  The 
output  can  be  varied  as  to  the  depth  and  detail  desired  f.  r 
a  given  input . 

The  inputs  used  for  the  ANMOL  program  are  also  those  used 
for  the  program  ANPLOT  which  will  provide  graphical  vice 
tabular  results.  In  addition  to  ANPLOT,  the  graphics  program 
EASYPLOT  was  used  to  generate  all  force  and  potential  plots. 

B.  HIT/NO  HIT  SCENARIO 

As  discussed  previously,  each  simulation  run  was 

conducted  with  a  total  of  300  impact  points  referenced  to  a 

pre-selected  point  on  the  target  surface.  For  this 

simulation,  only  half  of  the  available  positions  contain  a 

nitrogen  adatom.  Therefore,  there  is  an  equal  probability  of 

the  incident  ion  hitting  or  missing  an  adatom.  To 

incorporate  this  difference,  runs,  at  the  same  ion  energy, 

were  conducted  for  both  instances  where  the  ion  hit  and  did 

not  hit  a  nitrogen  adatom  with  the  results  being  averaged 

together.  This  results  in  essentially  a  600  point  run  for 

% 

each  energy  simulated.  The  primary  difference  of  the  two 
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runs  is  that  there  are  approxiraatly  two  to  four  times  the 
number  of  sputtered  nitrogen  atoms  for  the  hit  run  as  there 


are  for  the  no-hit  run,  depending  on  the  'incident  ion 
energy . 

All  table  listings  of  yield  and  cross  sections 
incorporate  this  hit/  no-hit  scenario. 


C.  NITROGEN  RESULTS 

1 .Simulated  cross  sections  and  yields 

The  value  of  the  sputtering  yield  of  nitrogen  can  be 
expressed  as  the  product  of  the  concentration  of  the 
nitrogen  adatoms  in  atoms  per  square  centimeter  (  0n  )  and 
the  sputtering  cross  section  in  square  centimeters  ( O' n  ) . 
The  results  are  in  units  of  nitrogen  atoms  per  ion.  Figure 
19  shows  the  experimental  results  of  Winters  for  the 
sputtering  cross  section  of  nitrogen  versus  the  ion  energy 
for  a  variety  of  substrates. 

The  simulated  concentration  of  the  nitrogen,  assuming  a 
coverage  of  0.4,  was  calculated  to  be  5.0E14  atoms /cm2 . 
Table  6  shows  the  results  of  the  nitrogen  cross  sections  for 
both  molybdenum  and  tungsten  substrates  with  the  adatoms 
located  0.05  LU  beneath  the  surface  plane. 

Since  the  adatom/substrate  potential  functions  change 
for  different  adatom  positions,  different  parameters  must  be 
calculated  for  each  assumed  position  of  the  adatom  to  ensure 
equilibrium.  New  potential  function  parameters  were  computed 


*.-4v - 


iy  .4'" ‘4  . 


■iW 


for  various  adatom  positions  above  the  plane  for  both  the 
molybdenum  and  tungsten  systems.  Table  7  shows  the  results 
for  the  selected  positions  ( in  parentheses  below  the 
substrate  type)  above  the  plane  for  each  system.  The  unit 
for  this  distance  is  the  LU. 


TABLE  6:  SIMULATION  NITROGEN  CROSS  SECTIONS 


CROSS 

.SECTION 

XIO-i s (cm2 

)  (BELOW 

PLANE) 

ION  ENERGY 

YIELD 

CROSS 

SECTION 

(KEV) 

Mo(001) 

W(001) 

Mo(001) 

W(001 ) 

1 

.  33 

N/A 

.  66 

N/A 

2 

.26 

.29 

.52 

.  57 

3 

.23 

N/A 

.46 

N/A 

TABLE  7:  SIMULATION  NITROGEN  CROSS  SECTIONS 

CROSS  SECTION  XlQ-is  (cm2 )  (ABOVE  PLANE) 


ION  ENERGY  YIELD  CROSS  SECTION 


(KEV) 

Mo(001) 

W(001) 

Mo(001) 

W(001) 

(  .245 

.38) 

(-487) 

(  .245 

.380) 

(  -487) 

1 

N/A 

.57 

.  567 

N/A 

1.14 

1.13 

2 

.301 

.  525 

.  450 

.6 

1.05 

.90 

3 

N/A 

.  55 

N/A 

N/A 

1.10 

N/A 

2.  Nitrogen  results  comparison  to  experimental  results 


Figure  20  shows  a  comparison  for  the  molybdenum 


substrate  with  adatom  positions  of  0.05  LU  below  the  surface 


for  the  low  point  and  0.38  LU  above  the  surface  for  the  high 
point. 


When  the  cross  section  and  yield  simulation  results  for 
the  molybdenum  target  with  the  below  plane  adatom  location 
are  compared  to  the  experimental  values  obtained  by  Winters, 


slight  positive  slope  Also,  the  simulation  cross  sr-r  1..t; 
values  are  below  the  experimentally  observed  values.  This 
indicates  that  the  adatom  probably  is  not  located  beneath 
the  surface  of  the  target.  The  observed  cross  section  for 
the  below  plane  nitrogen  from  tungsten  bombarded  at  2  keV 
was  obviously  low,  suggesting  that  the  adatom  was  not 
located  below  the  surface  for  that  system  either.  This 
observation  negated  the  need  for  further  runs  on  tungsten  at 
this  adatom  position. 

Figure  20  also  compares  the  computed  molybdenum  results 
for  the  two  previously  mentioned  above  plane  adatom 
locations  to  those  obtained  by  Winters.  With  respect  to  the 
0.38  LU  position,  the  curves  are  much  closer  in  shape,  and 
the  new  points  for  molybdenum  are  now  slightly  above  the 
experimental  values.  This  seems  to  indicate  that  the  assumed 
position  of  0.38  LU  above  the  plane  was  too  high  a  location 
for  the  adatoms. 

The  equidistant  positioning  of  the  adatom,  which  was  a 
change  of  0.135  down  LU  from  the  previously  examined  point, 
resulted  in  a  decrease  in  the  cross  section  of  40%.  This 
large  change  in  yield  for  a  small  change  in  distance 
supports  the  distance  dependency  postulated  by  Meyerhof f. 
These  results  also  indicate  that  the  position  of  the 
nitrogen  adatom,  for  a  molybdenum  target,  is  slightly  above 
the  eqidistant  position.  Figure  21  and  shows  a  comparison  of 


the  experimental  values  obtained  for  tungsten  substrate  and 


those  obtained  from  the  simulation  for  the  "above  plane' 


adatom  location  of  0.487  LU.  ’These  appear  to  be  in  fairly 


close  agreement.  However,  little  can  be  definitivly  stated 


since  only  two  points  were  evaluated  for  the  assumed 


position . 


3.  Simulated  nitrogen  energy  distribution 


Figures  22  through  25  show  the  energy  distribution  of  the 


sputtered  nitrogen  atoms  for  an  incident  energy  of  2  keV .  In 


each  case,  the  spectrum  was  fairly  broad.  In  both  extreme 


cases,  where  the  adatom  was  located  at  it’s  closest  and 


furthest  position  from  the  surface,  for  both  tungsten  and 


molybdenum,  a  spectrum  with  a  peak  of  approximaly  5  eV  is 


evident;  therefore,  the  choice*  of  the  substrate,  when  there 


is  a  large  mass  difference  between  the  substrate  and  the 


adatom,  appears  to  be  of  little  consequence  to  the  energy 


distribution  of  the  ejected  nitrogen.  Additionally,  the 


location  of  the  adatom  has  little  effect  on  it’s  sputtered 


energy.  These  results  conflict  with  those  obtained  by 


Meyerhof f . 


4. Simulated  nitrogen  ejection  angles 


Figure  26  is  a  spot  pattern  presentation  of  the  polar  and 


azmuthal  angle  dependence  of  the  sputtered  nitrogen  adatoms. 


The  azmuthal  angle  of  the  spot  is  read  normally  while  the 


polar  angle  of  a  spot  is  plotted  as  the  numerical  value  of 


the  tangent  of  the  polar  angle.  Therefore,  spots  further  out 


‘Tv 


-.c; 


correspond  to  adatoms  with  larger  angles  of  ejection  as 
measured  from  the  surface  normal.  Little,  if  any,  pattern 
development  can  be  seen  from  the  graph.  Table  8  lists,  in 
numerical  format,  typical  polar  and  azmuthal  angular  results 
for  sputtered  nitrogen  atoms  of  various  energies. 

Table  8*.  Nitrogen  atom  souttering  angles 
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D.  MOLYBDENUM  AND  TUNGSTEN  RESULTS 


When  the  incident  ion  impacts  a  target  with  a  chemisorbed 
gas;  target  atoms  as  well  as  the  adatoms  are  sputtered. 

Clean  surface  targets  of  both  molybdenum  and  tungsten 
were  bombarded  in  the  simulation  with  2  keV  argon  ions.  The 
number  of  sputtered  target  atoms  were  then  compared  to  the 
number  of  sputtered  target  atoms  when  the  target  had  been 
reacted  with  nitrogen.  Table  9  lists  these  values. 

TABLE  9:  SUBSTRATE  ATOM  EJECTIONS 


CLEAN 

ADATOM 

POSITION 

HIGH 

LOW 

907 

677 

686 

791 

590 

.  561 

A  comparison  of  the  number  of  sputtered  target  atoms  for 
the  low  and  high  nitrogen  positions  of  both  substrates  was 
done.  The  comparison  shows  that  although  the  nitrogen  acts 
as  an  energy  absorption  mechanism.  it's  location  is  of 
little  consequence  to  the  ejection  mechanism  of  the 
substrate  material. 

1 .  Simulation  target  yield 

Figures  27  and  28  show  the  target  atom  yield  per  impact 
point  for  a  2  keV  incident  ion.  The  bottom  left  corner  of 
the  diagram  is  the  selected  reference  point  for  the  ion 
impact.  Throughout  the  graph,  the  darker  the  point,  the 
higher  the  yield  for  that  particular  point.  The  lighter 
areas  of  the  graphs  are  indicative  of  the  channelling  effect 
for  the  BCC  structure. 
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a.  Simulation  Target  Yield  Comparison  to  Experiment 
In  his  experiments,  Winters  also  calculates  the 
ratios  of  the  cross  sections  of  the  chemisorbed  gas  to  the 
substrate  material.  Table  10  lists  the  results  of  the 
simulation  as  compared  to  the  experiment.  Winters  used 
polycrystalline  tungsten  and  molybdenum  targets  in  his 
experiments  while  the  calculations  of  this  thesis  used  clean 
target  (no  adatoms)  cross  sections.  The  simulation  tungsten 
and  molybdenum  concentrations  were  calculated  to  be  1.19E15 
and  1.2E15  respectivly. 


TABLE  10:  CROSS  SECTION  RATIO’S  WITH  NITROGEN  IN  EQUIDISTANT 
POSITION 


EXPERIMENT 

ION  ENERGY  <TN(W)  C'.Nl.MQJL- 

CTW(POLY)  O' MO  ( POLY ) 
2  .77  .68 


SIMULATION 
g^N(W)  CT  N  ( MO ) 

<TU(  CLEAN)  O' MO  (CLEAN) 
.41  .24 


f 


Since  the  simulation  clean  targets  were  bombarded  with  2  keV 
ions  only,  comparisons  must  be  restricted  to  that  energy. 
The  simulation  cross  section  value  for  nitrogen  from  the 
molybdenum  target  was  calculated  assuming  the  0.245  position 
for  the  adatom.  Table  10  shows  a  smaller  value  than  the 
experimental  ratio  although  the  nitrogen  cross  section  ratio 
for  experiment  and  simulation  are  both  1.33.  These  results, 
along  with  the  previously  mentioned  significant 
yield/distance  relationship  indicate  that  the  assumed 


positions  are  close  to  the  actual  adatom  locations. 


2 .Simulation  probabilty  of  ejection 

Figures  29  and  30  are  representitive  of  the  probability 
of  ejection  for  each  of  the  surface  atoms,  in  this  case  on 
the  (001)  face.  These  figures  indicate  how  well  the  cascade 
sequence  was  contained  by  the  target.  The  graphs  show  that 
the  size  of  the  target  was  generally  acceptable  for 
molybdenum  but  was  slightly  small. for  tungsten.  This  does 
not  affect  the  sputtering  results  for  either  the  substrate 
or  the  adatom  significantly. 

3.  Simulation  substrate  ejection  spot  patterns 

Figures  31  and  32  show  typical  sputtered  target  atom 

spot  patterns.  As  with  the  nitrogen  graph,  these  give 

graphical  representations  of  the  angles  of  ejection  of  the 
target  atoms.  The  crystalline  nature  of  the  targets  is 
readily  apparent  in  these  graphs  because  of  the 
concentration  of  sputtered  atoms  along  the  45® , 135° , 225° , 
and  315°  arms.  For  a  BCC  structure,  these  directions 

correspond  to  the  nearest  neighbor  directions 
4 . Simulation  ejected  substrate  energy  distribution 

Figures  33  through  38  show  typical  energy  distribution 
histograms  of  the  sputtered  target  atoms.  There  were  no 
apparent  differences  between  the  specific  substrate  results. 
The  molybdenum  results  show  a  maximum  at  approximatly  5eV 

while  the  maximum  for  tungsten  is  slightly  higher.  These 
values  conflict  with  the  theoretical  values  which  should  be 
one  half  the  cohesive  energy;  specifically,  3.41  eV  for 


molybdenum  and  4.45  eV  for  tungsten.  When  plotted 
logarithimicly ,  both  elements-  approximatly  show  the  expected 
E- 2  dependency  at  higher  energies  (Fig. 39  through  43).  Table 
11  shows  the  slope  relationship  numerically  with  the  values 
having  been  obtained  from  the  figures  mentioned. 

TABLE  11:  SUBSTRATE  EJECTED  ATOM  ENERGY  DISTRIBUTION 
SLOPES 


ION  ENERGY 
(KEV) 

SUESTRATE 

ENERGY 

DISTRIBUTION  SLOPE 

3.0 

MO 

-1 . 8 

1.0 

W 

-1.93 

2.0 

w 

-1 . 94 

2.0 

W  (clean) 

-1 . 99 

2.0 

MO  (clean) 

-2.13 

E.  DIMER  RESULTS 

When  a  compound  structure  is  bombarded,  there  may  be  an 
expectation  that  groups  of  two  or  more  atoms  would  be 
sputtered  off  together.  The  results  of  this  simulation  do 
not  show  this  to  be  the  case.  In  many  instances  however, 
after  ejection,  two  atoms,  either  substrate  -  substrate  or 
substrate  -  adatom  combine  over  the  target.  Individual 
ejection  times  of  combining  atoms  showed  differences  of  as 
much  as  200  fsec,  indicating  that  combination  occured  at  a 
relativly  large  distance  from  the  target.  However,  there 
were  some  instances  when  the  times  of  ejection  were 
separated  by  as  little  as  10  fsec  for  two  atoms  whose 
pre-ejection  locations  were  within  1  LU  thus  indicating 
dimer  formation  almost  immediatly  after  ejection. 
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7.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  results  of  the  simulations,  when  compared  to 
experimental  results,  indicate  that,  for  nitrogen 
chemisorbded  on  the  (001)  face  of  molybdenum  and  tungsten 
targets,  the  location  of  the  adatom  is  above  the  surface 
plane  for  both  targets. 

Nitrogen  yield  comparisons  between  the  simulation  and 
Winters  experimental  results  indicate  that  the  adatom 
position  for  the  molybdenum  target  is  above  the  equdistant 
position  of  0.245  LU  and  below  the  0.380  LU  position. 
Furthermore,  the  simulation  yield  comparisons  between  the 
two  positions  show  that  there  is  a  very  strong  dependency  of 
the  yield  on  the  adatom  distance  from  the  target  surface 
plane . 

The  comparison  of  the  simulation  nitrogen  yield  to  the 
experimental  values  for  the  tungsten  target  indicate  that 
the  adatom  position  is  near  the  equidistant  position  of 
0.487  LU.  However,  since  only  two  ion  energy  values  were 
analyzed  for  this  target,  further  investigation  of  the 
tungsten  target  should  be  conducted  at  different  ion 
energies  to  complete  the  argon  simulation  results. 

Additionally,  further  simulation  investigation  into  both 
molybdenum  and  tungsten  should  be  conducted  using  different 
bombarding  ions  in  the  same  energy  range  as  those  used  in 
this  thesis.  These  additional  simulations  should  be  analyzed 


WWW 


in  conjunction  with  the  results  obtained  in  this  thesis  to 
isolate  the  exact  adatom  locations  for  the  targets  in 
question . 
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Small  circle  represents  nitrogen  atoms. 

Fig.  1.  Location  of  nitrogen  atoms  on  target  substrates 


The  proposed  contracted-domain  structure 
for  the  0.4  monolayer  on  a  Mo(OUJ)  or  a 
W(001)  surface:  Large  hatched  circles,  top- 
layer  substrate  atoms,  small  filled  circles 
N  atoms. 

a)  Plan  view,  illustrating  domain  and  boundary 
structure. 

b)  Cross  section  through  the  domain  along  AA: 
Fig.  2.  Island  formation  on  a  target  surface 
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MO-N  FORCE  (LOW  POSITION) 


Fig.  5.  Molybdenum-Ni trogen  force  function,  adatom  location 
-0.05  LU 
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Fig.  7.  Molybdenum-Nirtogen  potential  function 
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W-N  POTENTIAL 


■  0.5  1.0  1.5  2.0 


SEPARATION  (LU) 
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Fig.  12.  Molybdenum-Molybdenum  force  function 
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SEPARATION  (LU) 

Fig.  14.  Tungsten-Tungsten  force  function 
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N-N  POTENTIAL 


15.  Nitrogen-Nitrogen  potential  function 
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10  12  14  16  18 


a)  A-top  position 


b)  2  fold  position 


Fig.  18.  Nitrogen  adatoms  in  various  positions 
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Fig.  20.  Simulati  on  nitrogen  sputtering  cross  sections 
for  molybdenum  target 
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W  SPUTTERING  CROSS  SECTIONS 


ION  ENERGY  (KEV) 


Fig.  21.  Simulation  nitrogen  sputtering  cross  section 
for  tungsten  target 
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FRACTION 

0.00  0.05  0.10  0.15  0.20 


2.000  W(001)  +60N/A<001>  COMBINED  (N)(LOW) 

EJECTED  ATOM  ENERGY  DISTRIBUTION 


Fig.  22.  Ejected  nitrogen  atom  energy  distribution  iron 

a  tungsten  target  assuming  Argon  energy  of  2  keV 
and  adatom  location  of  -.05  LU 
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2.000  MO (001)  +60N/A<001>  COMBINED  (N)(LOW) 

EJECTED  ATOM  ENERGY  DISTRIBUTION 


ENERGY  (EV) 


Fig.  23.  Ejected  Nitrogen  energy  distribution  from  a 

molybdenum  target  assuming  argon  energy  2  keV 
and  adatom  location  -.05  LU 


2.000  M0(001)  +60N/A<001>  COMBINED(N)(HIGH) 

EJECTED  ATOM  ENERGY  DISTRIBUTION 


ENERGY  (EV) 

Fig.  24.  Ejected  nitrogen  energy  distribution  from 

molybdenum  target  assuming  argon  energy  of  2  keV 
and  adatom  position  . 387  LU 
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Fig.  25.  Ejected  nitrogen  atom  energy  distribution  from  a 
simulated  tungsten  target  assuming  argon  energy 
of  2  keV  and  adatom  position  .487  LU 
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2.000  M0(001)  +60N/A<001>  NITROGEN  (NH)(HIGH) 
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Fig.  26.  Simulation  ejected  nitrogen  spot  pattern  from  a 
molybdenum  target  assuming  argon  energy  of  2  keV 
and  adatom  location  . 387  LU 
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2.000  W(001)  +60N/A<001>  (23X8X23)  TUNGSTEN  (NH) 
ATOM  YIELD  PER  IMPACT  POINT 
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Fig.  29.  Simulation  atom  ejection  probability  from  a 
PCC  (001)  molybdenum  target  assuming  a  2keV 
argon  ion 
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2.000  W(001)  +60N/A<001>  (23X8X23)  TUNGSTEN  (NH) 

ATOM  EJECTION  PROBABILITY 

BCC  (001) 


Fig.  30.  Simulation  atom  ejection  probability  from  a 
BCC  (001)  tungsten  target  assuming  an  argon 
ion  energy  of  2  keV 
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2.000  W(001)  +60N/A<001>  TUNGSTEN  (NH)(LO\V) 


Fig.  32.  Simulation  tungsten  ejection  spot  pattern  assuming 
2  keV  argon  ion 
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Fig.  33.  Simulated  ejected  atom  energy  distribution  from 
molybdenum  target  with  nitrogen  adatoms  assuming 
a  3  keV  argon  ion 
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Fig.  34.  Simulated  ejected  atom  energy  distribution  from 
tungsten  target  with  nitrogen  adatoms  assuming 
a  2  keV  argon  ion 
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Fig.  35.  Simulation  ejected  atom  energy  distribution  from 
tungsten  target  with  nitrogen  adatoms  assuming  a 
1  keV  argon  ion 
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Fig.  36.  Simulation  ejected  atom  energy  distribution  from 
molybdenum  target  with  nitrogen  adatoms  assuming 
a  2  keV  argon  ion 
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Fig.  37.  Simulated  ejected  atom  energy  distribution  from 
clean  molybdenum  target  assuming  a  2  keV  argon 
ion 
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Fig.  38.  Simulated  ejected  atom  energy  distribution  fro 
clean  tungsten  target  assuming  a  2  keV  argon  i 
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Fig.  39.  Log-Log  ejected  atom  energy  distribution  from  a 
molybdenum  target  with  nitrogen  adatoms  assuming 
3  keV  argon  ion 
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Fig.  40.  Log-Leg  ejected  atom  energy  distribution  from 
tungsten  target  with  nitrogen  adatoms  assuming 
1  keV  argon  ion 
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Fig.  41.  Log-Log  ejected  atom  energy  distribution  from 
tungsten  target  with  nitrogen  adatoms  assuming 
2  keV  argon  ion 
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Fig.  42.  Log-Log  ejected  atom  energy  distribution  from 
clean  tungsten  assuming  2keV  argon  ion 
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Fig.  43.  Log-Log  ejected  atom  energy  distribution  from 
clean  molybdenum  assuming  2  keV  argon  ion 
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